Spin-triplet superconductivity in repulsive Hubbard models with disconnected Fermi 
surfaces: a case study on triangular and honeycomb lattices 
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We propose that spin-fluctuation-mediated spin-triplet superconductivity may be realized in 
repulsive Hubbard models with disconnected Fermi surfaces. The idea is confirmed for Hubbard 
models on triangular (dilute band filling) and honeycomb (near half-filling) lattices using fluctuation 
exchange approximation, where triplet pairing order parameter with f-wave symmetry is obtained. 
Possible relevance to real superconductors is suggested. 

PACS numbers: 74.20.Mn, 71.10.Fd 



A fascination toward spin-triplet superconductivity 
has a long history, but recent experimental suggestions 
for triplet pairing in a heavy fermion system UPt3, JlJ or- 
ganic conductors (TMTSF) 2 X (X=C10 4 [§,PF 6 §]), and 
a ruthenate Sr 2 Ru04, |4j|5|] have renewed our interests 
in mechanisms of triplet superconductivity. In particu- 
lar, it is fairly intriguing to investigate whether electron- 
electron repulsive interactions can lead to triplet super- 
conductivity. Ferromagnetic-spin-fluctuation mechanism 
has been proposed from the early days, but to our knowl- 
edge, realization of triplet superconductivity (at sizable 
temperatures) has yet to be established theoretically for 
repulsive electron models with renormalization effects of 
the quasiparticles taken into account. The lifetime of the 
quasiparticles is important since this is a factor dominat- 
ing T c . 

Recently, the present authors with Aoki have investi- 
gated the possibility of triplet pairing in the Hubbard 
model for a variety of lattice structures and band fillings 
using fluctuation exchange (FLEX) approximation. || A 
naive expectation is that triplet superconductivity may 
be realized when the band is away from half-filled and 
the density of states (DOS) at the Fermi level is large, 
since ferromagnetic fluctuations become strong in such 
a situation. In ref. ||, however, it has turned out that 
the transition temperature (T c ) of triplet superconduc- 
tivity, if any, is too low to be detected as far as the cases 
surveyed there are concerned. A typical case is a square 
lattice with appreciable next nearest neighbor hoppings 
and dilute band fillings. First let us briefly review this 
situation as a reference for the results presented later. 



We 



consider 
t-c 1 ' 



the Hubbard model, Tt = 
12(i,j)<j=7,itij c i<T C j<? + Uj2i n iT n il, on a square lattice 
shown in Fig.[l], where t(= 1) is the nearest and t[(= t' 2 
here) is the next nearest neighbor hopping. In the FLEX 
calculation, raj (i) Dyson's equation is solved to ob- 
tain the renormalized Green's function G(k), where k 
k. ie n ) denotes the 2D wave- vectors and the Matsubara 
frequencies, (ii) the effective electron-electron interaction 
V^\q) is calculated by collecting RPA-type bubbles and 



ladder diagrams consisting of the renormalized Green's 
function, namely, by summing up powers of the irre- 
ducible susceptibility Xirr(<z) = — jj J2k G(k + q)G(k) 
(iV:number of fc-point meshes), (iii) the self energy is 
obtained as E(fc) = i J2 q G ( k ~ <?)^ (1) (<?), which is 
substituted into Dyson's equation in (i), and the self- 
consistent loops are repeated until convergence is at- 
tained. Throughout the study, we take 64 x 64 k- 
point meshes and the Matsubara frequencies e n from 
~(2N C - 1)ttT to (2N C ~ 1)ttT with N c up to 16384 in 
order to ensure convergence at low temperatures. 

To obtain T c , we solve the eigenvalue (Eliashberg) 
equation for the superconducting order parameter </>(fc), 



X ^ k ) = ~7j E ^{k')\G{k')\ 2 V^{k k% 



N 



(1) 



where the pairing interaction , which mediates pair 
scattering from (k, — k) to (k', — k')[= (k + q, k — q)j, 
is given as 



1 - UXirM 



for singlet pairing, and 



u 2 xM 



UXirr{q) 



U 2 XirM 



1 + UXirM 



(2) 



(3) 



for triplet pairing. In either case, the first (second) term 
represents the contribution from spin (charge) fluctua- 
tions. T c is the temperature at which the maximum 
eigenvalue A reaches unity. We denote the eigenvalue 
and the order parameter for triplet (singlet) pairing as 
A t (A s ) and <f> t (cj> s ), respectively. 

In ref. [§, t', U, and the band filling n [|| were varied 
in search of triplet superconductivity, but At remained 
below ~ 0.2 in the tractable temperature range as typ- 
ically displayed in Fig.|l| (dash-dotted line). The main 
reason why triplet pairing instability is weak is because 

(2) (2) 

\V t I is only one third of \Vs | when spin fluctuation is 
dominant as can be seen from eqs.(^) and (||). 
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FIG. 1. The left panel shows the hopping integrals for 
square (t[ — t' 2 ) or triangular (t' 2 = 0) lattice. In the right 
panel, At is plotted as a function of temperature for the Hub- 
bard model on a square lattice with t'i — t' 2 — 0.4, (7 = 6, 
and n = 0.3 (dash-dotted line), or on an isotropic triangular 
lattice with t[ = 1, U = 8, and n = 0.15 (solid line). In 
the latter spline extrapolation to lower temperatures 

is also plotted (dashed line). 

In this Letter, we propose that the above difficulty 
for spin-fLuctuation-mediated triplet pairing in the Hub- 
bard model can be overcome under certain conditions. 
Let us first present our idea. We consider a situation 
(see Fig.||) where (i) the Fermi surface (FS) is discon- 
nected (preferably well separated) into two pieces which 
are located point symmetrically about k = 0, and (ii) 
the spin structure is pronounced around a wave vector 
Q in such a way that two electrons with zero total mo- 
mentum can be scattered within each piece of the FS 
(this process will be called intra-FS pair scattering here- 
after) by exchanging spin fluctuations having momen- 
tum ~ Q. Now, in order to have large A, the quantity 

-E k , k ^F5^ (2) (k-k')^(k)0(k')]/E ke i,5^( k )] ( thc 

numerator being ~ — X^ke-FS ^ 2 HQ)^(k)<Hk + Q)) has 
to be positive and large. Then, pair scatterings from 
(k, — k) to (k + Q, — k — Q) for singlet pairing, mediated 

f 2) 

by repulsive V s (Q), have to accompany a sign change in 
the order parameter </5 s (k) (Fig.||(a)). Hence the nodes 
of <j) s (k) must intersect the FS. For triplet pairing, by 
contrast, pairs can be scattered within a region having 

( 2 ) 

the same sign in 0t(k) because V t (Q) is attractive. In 
this case, since the gap nodes (which exists due to triplet 
pairing symmetry cf>(k) — — <fi(— k)) do not intersect the 
FS (Fig.||(b)), the entire FS can be exploited for pairing, 




FIG. 2. Basic idea of the present mechanism is presented. 
Intra-FS pair scattering (dashed curves) is mediated by 
V (2) (Q) within each pieces of the disconnected FS (closed 
solid curves). The dashed straight lines represent nodes of 
the order parameter for singlet (a) and triplet (b) pairing. 




FIG. 3. \G{K,mk B T)\ 2 (a), x(k,0) (b), and <f> t {V.,mk B T) 
(c) are plotted for the Hubbard model on a triangular lattice 
with t[ = 1, U = 8, n = 0.15, and T = 0.01. The right panels 
are contour-plots of the left panels. 

so that triplet pairing may be enhanced. Quite recently, 
a related proposal has been raised by Kohmoto and Sato 
for systems with both phonons and spin fluctuations 
present, B as discussed later. 

The above conditions are not satisfied for the t-t' 
square lattice because it has a connected FS. As an ex- 
ample of a system in which the above conditions are 
indeed satisfied, we consider the Hubbard model on 
an isotropic triangular lattice with dilute band fillings. 
The band dispersion for U — is given by £A(k) = 
2[cosfca; + cos k y + cos(k x + k y )\ when we represent an 
isotropic triangular lattice by setting t = t[ = 1 M] and 

= in Fig^. Superconductivity on an isotropic trian- 
gular lattice has been studied by several authors, |ll] M 
but their interest was mainly focused on n ~ 1. In ref 
||, possibility of triplet superconductivity was studied 
at quarter filling (n = 0.5), where ferromagnetic fluctua- 
tions become strong because the Fermi level lies right at 
the position where the DOS diverges. However, At was 
again found to be small, which, in the present context, 
is because the FS is connected. If we set n < 0.5, on 
the other hand, the FS is disconnected into two pieces, 
which are centered respectively at k = (2ir/3,2w/3) and 
k = (47r/3, 4-7r/3). Here we take n = 0.15, where the two 
pieces of the FS are well separated. In Fig.||, we plot the 
FLEX result for |G(k, inksT)] 2 (a) and the spin suscep- 



tibility x (k,0) = Xhr(k,0)/[1 - C/x irr (k,0)] (b) against 
k for U = 8 and T = 0.01. The FS as identified from 
the ridge in |G(k, i7rfcsT)| 2 is indeed disconnected into 
two pieces. %(k, 0) is sharply peaked at k = as seen in 
Fig.||(b), indicating ferromagnetic fluctuations. [p| This 
is partially because the FS is small, but it is also because 
the Fermi level for n = 0.15 is still not too far away from 
the peak position of the DOS. In this case, A s is shown 

(21 

to be small, which is because V s (Q) can only mediate 
pair scatterings in the vicinity of the nodes when Q ~ 0. 

If we turn to triplet pairing, the order parameter 
<£ t (k, wrfc B T), plotted against k for T = 0.01 in Fig.|(c), 
has f-wave (f a ,3_ 3;1 , J( 2-wave in the notation of the Cq sym- 
metry group) symmetry with three sets of nodal lines 
(k x = 0(mod27r), k y = 0, and k x + k y = 0). Compar- 
ing Figs.||(a) and (c), we can see that these nodes do 
not intersect the FS. Accordingly, At (Fig.|l|,solid line) 
is strongly enhanced compared to the case for the t-t' 
square lattice. A spline extrapolation to low tempera- 
tures suggests a possible low but finite T c . 

As another example, we next propose that the Hub- 
bard model on a honeycomb lattice (Fig.[|) should also 
be interesting. Since there are two sites (A and B) in 
a unit cell, this is a two-band system. The noninter- 
acting band dispersion £i lc (k) = ±^/eA(k) + 3 has two 
pairs of vertex-sharing cones at k = (27r/3, 2ir/3) and 
k = (47r/3, 47r/3), so again the FS becomes disconnected, 
this time for fillings close to n = 1. 

In the multiband version of FLEX, HQ the quan- 
tities G, x, X, and 4> have 2x2 matrix forms, e.g., 
G Q/ 3(k, iu n ), where a, (3 denote A or B sites. The band 
representation of the Green's function and the order pa- 
rameters is obtained by using the relation between the 
annihilation operators of upper (u) and lower (I) band 
electrons (c£, c^.) and those of A and B site electrons (c£, 
c^). As for x, we diagonalize the 2x2 matrix Xa/3 to ob- 
tain x± = (xaa+Xbb)/2±\/[(xaa - Xbb)/2] 2 + |xab| 2 - 

In FigJ| we plot \G l (k,mk B T)\ 2 (a) and x+(k,0) (b) 
for n = 0.95, U = 8 and T = 0.01. Since xab(0,0) is 
found to be negative, the peak around k = in X+(k, 0) 
is an indication of antiferromagnetic fluctuations, as ex- 
pected for a nearly half-filled bipartite lattice system. 
Note that X+(k, 0) has a broad structure compared to 
the case for the triangular lattice (Fig|3|(b)). 




FIG. 4. For the honeycomb lattice shown in the left panel, 
we employ the topologically equivalent structure shown in the 
right. 
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FIG. 5. \G l (k,i7vk B T)\ 2 (a), x +(k,0) (b), $(k, ink B T) (c), 



and 03 (k, i^ksT) (d) are plotted for the Hubbard model on 
a honeycomb lattice with U = 8, n = 0.95, Aab = and 
T = 0.01. 

If we turn to X s and At as functions of T in Fig.^|(a), 
At is again large, but this time A s is in fact larger. We 
can trace this to the broad spin structure, for which spin 
fluctuations with relatively large momentum can be ex- 
changed to mediate singlet pairing at wave vectors away 
from of the nodes. Nevertheless, we can still observe that 

(21 (2) 

At is enhanced above A s /3 (recall that \V t \ ~ \Vs |/3), 
which should be due to the fact that the nodes in <p l s 
intersect the FS, while those in <j)\ do not as seen by 
comparing Figs.||(a) and (c)/(d). 

We have found that |<^ab| > (<) |</"aa| f° r singlet 
(triplet) pairing, meaning that singlet (triplet) pairing 
mainly takes place on different (same) sublattices. This 
is in fact consistent with the antiferromagnetic alignment 
of the spins. Then, we can intuitively expect that triplet 
can dominate over singlet if we introduce a level offset, 
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FIG. 6. At (solid line) and A s (dash-dotted line) are plot- 
ted as functions of temperature for the Hubbard model on 
a honeycomb lattice with U = 8, n — 0.95, and Aab = 
(a) or Aab = 5 (b). In (b), spline extrapolations to lower 
temperatures are indicated by dashed lines. 

Aab, between A and B sites. In Fig.^(b), we plot A s 
and A t for the same parameter values as in Fig.^(a), ex- 
cept for a finite Aab = 5. Triplet pairing indeed domi- 
nates over singlet pairing at low temperatures, and here 
again a possible finite T c for triplet superconductivity is 
suggested. The intuitive picture can be paraphrased in 
the momentum space that the peak structure of x ( n °t 
shown) becomes sharper when Aab ^ 0. 

Finally, let us make some remarks concerning possi- 
ble relevance to real superconductors. As for UPt3, if 
we examine the FS calculated by FLAPW method, 
we notice that there are two disconnected pockets (band 
37 in ref. g§), @ which, in our view, is favorable for 
triplet pairing. It would be an interesting future prob- 
lem to investigate in detail the applicability of the present 
mechanism. 

Disconnected FS can arise similarly in graphite interca- 
lation compounds (GIC), except that the FS is cylindri- 
cal (quasi 2D), [^(j This is because graphite is a system 
where honeycomb sheets of carbon atoms are stacked. 
Although spin fluctuations in GIC may not be strong 
enough to induce superconductivity purely electronically, 
the disconnectivity of the FS itself should be favorable for 
triplet pairing, so a cooperation between certain phonon 
modes and (weak) spin fluctuations might lead to triplet 
superconductivity (even in the absence of Aab consid- 
ered above) . Namely, if attractive intra-FS pair scatter- 
ings mediated by phonons are present, antiferromagnetic 
spin fluctuations as considered here would work construc- 
tively with phonons to enhance intra-FS pair scatterings 
for triplet pairing, while the converse is true for singlet 
pairing. Experimentally, although triplet pairing has not 
been claimed to our knowledge, a large value of H C 2 (ex- 
trapolated to T = 0) observed in C4KTI1.5 ]2l]] is in fact 
reminiscent of a large H c2 in (TMTSF) 2 X. ||| 

As for (TMTSF) 2 X and Sr 2 Ru0 4 , Kohmoto and Sato 
have recently proposed that disconnectivity (quasi-one 
dimensionality) of the FS, along with the presence of spin 
fluctuations originating from the nesting of the FS, plays 
an essential role in stabilizing phonon-mediated triplet p- 
wave pairing. M Our study is related to this proposal 



in that disconnectivity is important, but in these sys- 
tems, as seen in Kohmoto and Sato's argument, the dom- 
inant spin fluctuations have wave- vectors that bridge the 
two pieces of the FS, so that they mediate inter-FS pair 
scatterings p3] rather than intra-FS ones. Thus, our 
purely spin-fluctuation-mediated pairing mechanism does 
not directly apply to these materials, although we do be- 
lieve that the disconnectivity of the FS may be playing a 
certain role in realizing triplet superconductivity. 
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